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EXPERIMENTAL DETERMINATION OF AIRPLANE 

MASS AND INERTIAL CHARACTERISTICS 

Chester H . Wolowicz and Roxanah B. Yancey 
Flight Research Center 


INTRODUCTION 


Mass and inertial characteristics— the location of the center of gravity , the 
moments of inertia , and the inclination of the principal axis— must be known to pre- 
dict the performance and limiting flight conditions of new airplane configurations . 

The accuracy of such predictions is particularly important for unconventional ve- 
hicles like the lifting bodies . Furthermore , after a flight test program begins , the 
accuracy of the predicted mass and inertial characteristics affects the precision 
with which the stability and control derivatives can be determined from the flight 
data. 

Estimates of airplane mass and inertial characteristics are made and refined 
during the design process . To verify the calculated characteristics , the mass dis- 
tribution characteristics of a new configuration are frequently determined experi- 
mentally . However , in any experimental investigation there are limitations on the 
experimental setups and precautions to be observed . In these respects , the experi- 
mental techniques have continuously improved . 

At first, experimental mass characteristics were determined from overhead sus- 
pensions (refs . 1 to 4) . The vertical center-of-gravity position was obtained by 
tilting the airplane through large pitch angles and sighting the point of intersection 
of the line of suspension in the XZ -plane. The pitching and rolling moments of 
inertia were obtained by suspending the airplane as a compound pendulum . The 
yawing moment of inertia was obtained with the airplane suspended as a bifilar 
pendulum . The inclination of the principal axis was obtained from the experimen- 
tally determined ellipsoid of the inertias . However , structural , handling , and accu- 
racy problems make it impractical to apply these methods to airplanes that weigh 
more than approximately 45,000 newtons. 

Improved and simpler techniques were reported in reference 5 , The pitching 
and rolling moments of inertia were obtained by the compound-pendulum method , 
but the airplane was supported on and pivoted about knife edges , and the restoring 
moments were provided by springs . The yawing moment of inertia and the inclina- 
tion of the principal axis were determined by suspending the airplane support system 
from a vertical torsion tube in a way that permitted freedom in pitch but restrained 



the airplane's rolling motion by the bending of the tube. Reference 5 is the first 
reported instance of finding the roll-to-yaw ratio of the response to yaw excitation 
as a function of pitch attitude to obtain the inclination of the principal axis . Evi- 
dently , there were practical problems that precluded further use of the torsion 
tube, such as difficulty in providing rigid restraint at the fixed end of the tube. 

In current practices of determining mass and inertial characteristics , the ver- 
tical center-of- gravity position is obtained by tilting the airplane through a range 
of pitch attitudes while it is mounted on weighing scales located at jack points or 
the landing gear . Reference 6 considers this method to be accurate to within 5 per- 
cent . An overhead suspension method that provides more precise results is in use 
at the NASA Flight Research Center, To determine pitching and rolling moments of 
inertia, the airplane is supported on and pivoted about knife edges. To determine 
yawing-moment inertias and the inclination of the principal axis , a single-point sus- 
pension double-pendulum technique is used . The elementary principles of this 
technique are described in reference 7 . Reference 8 expands on the discussion in 
reference 7 and provides a comprehensive theoretical treatment to insure proper rig 
design for the condition where the restoring-moment springs are in planes parallel 
to the YZ-plane . At the NASA Flight Research Center the same technique is used, 
but the springs are in planes parallel to the XZ -plane. 

Although the literature has been concerned with improved techniques for the 
experimental determination of airplane mass and inertial characteristics , little 
mention has been made of the effects of amplitudes of oscillation on the degree of 
nonlinearity in the gravity-foree and spring moments in pitching- and rolling- 
moment-of-inertia determination . 

This report reviews current practices for experimentally determining the air- 
plane center of gravity, moments of inertia, and products of inertia. The techniques 
discussed are applicable to bodies other than aircraft. The method used at the NASA 
Flight Research Center to obtain the vertical center-of-gravity position accurately is 
described , along with the more conventional method which involves tilting the air- 
plane at jack points . The report provides information to facilitate the assessment of 
proposed experimental setups for pitching- and rolling-mom ent-of-inertia tests in 
terms of the amplitude of oscillations permissible for the accuracy desired. Finally, 
the NASA Flight Research Center’s application of the single-point suspension double- 
pendulum method to determine yawing moment of inertia and product of inertia is 
described, and the interaction of the modes of oscillation and the requirements to 
suppress or minimize the undesired modes are reviewed . 


SYMBOLS 


Physical quantities in this report are given in the International System of 
Units (SI) . Details concerning the use of SI are given in reference 9 . 


A = 



N 


2 



a 


B = ± 



length of the line of action of the spring , in yawing-moment- 
of-inertia tests , from the attachment point at the aircraft 
to the tieback stanchion (fig. 11 ) , m 


= *"('3)' (♦ ®sp ' ‘Pi ®sp)’ 



d 


horizontal distance between the front and rear loading points 
used in obtaining the horizontal center-of-gravity position 
by the methods shown in figures 1 (a) , 1 (b) , and 2 , m 


E 


s 



g 


horizontal distance between the center-of-gravity datum 
and the rear loading point used in obtaining center-of- 
gravity positions by the method shown in figure 2 , m 

ratio of the nonlinear increment in the change in spring 
moment to the linear change for the aircraft oscillating 
about a pivot 

ratio of the nonlinear increment in the change in weight 
moment to the linear change for the inverted-pendulum 
concept of the aircraft oscillating about a pivot 

acceleration due to gravity , m/sec^ 




moments of inertia of the aircraft about the X- , Y- , and 
Z-body axes, respectively, kg-m^ 



moment of inertia of the aircraft about its Y-axis passing 
through its center of gravity , kg-m^ 




moments of inertia of a test-setup component (such as 
cradle) about the Y- and Z-axes passing through the 
component's center of gravity, respectively, kg-m^ 



moments of inertia of a test-setup component 
(such as cradle) about the X- and Y-pivots of 
oscillation, respectively, kg-m^ 


Iv-7 product of inertia relative to aircraft X- and Z-body axes , 

1 2 
kg-m 

z product of inertia of an aircraft relative to the X- and Z- 

, s s axes passing through the line of suspension , kg-m^ 
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I- 5 J. rj with test-setup components included , kg-m^ 


product of inertia of a test-setup component relative to the 
X- and Z-axes passing through the component's center 
of gravity , kg-m^ 

product of inertia of a test-setup component relative to the 
X- and Z-axes passing through the line of suspension , 
kg-m^ 



change in pitching moment of inertia due to apparent 
additional air mass effects , kg-m^ 




yawing moment of inertia of an aircraft about the line of 
suspension, kg-m^ 

change in yawing moment of inertia of an aircraft about 
the line of suspension due to apparent additional air 
mass effects , kg-m^ 

yawing moment of inertia of an aircraft in the test setup 
about the line of suspension , kg-m^ 

linear spring constant of an individual spring , N/m 


k^,k 2 ,k„,k. spring constants of springs used in single-point suspen- 

sion test setup (fig . 17) , N/m 

k^ torsional spring constant, N-m/rad 

k = ir ' 
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t sp 


N-m/rad 




L ,N ,Y 
s s s 


2 2 
a/c’^W 


sinusoidal forcing rolling-moment input , N-m, and side 
force input , N, respectively 

spring rolling moment, N-m, yawing moment , N-m, and 
side force , N , respectively 

distance from the center of gravity of airplane and com- 
ponent weight , respectively , to the pivot of oscillation 
normal to pivot axis , m 

length of suspension beam (used in appendix B) , m 


4 



£ £ £ 
0’ 1’ 2 


£. = £^ cos 6- , m 

ih i i 

£'' £'' £'' 

0’ 1’ 2’ 3 


AM 


AM 


W 






w , 

m = g- , kg 

""a/c ^ 

"^a/c = — ’ 

W. 


t> 1 


P 1,P 2,P 3 
n n n 


P 1,P 2 

V V 


distance between reference rivets , defining reference 
longitudinal body axis (fig . 3) , m 

dimensions used in pitch and roll configurations , defined 
in table 1 and figure 6 , m 


dimensions used in single-point suspension method for 
determining I„ , I„_ , and e , defined in figures 11(a) 
and 13, m ^ 

change in spring moment, N-m 

change in weight moment due to perturbation of inverted 
pendulum action of weight , N-m 

AM^ when pendulum arm is horizontal , N-m 


period of oscillation , see 

spring configurations in pitch tests with springs normal to 
moment arm extending from spring attachment point on 
airplane to pivot of oscillation (fig . 6) 

spring configurations in pitch tests with springs vertical 


P.P 

|pM 
("^|p|)rm 
IPl 


|r| 




Pi = 


roll rate, rad/sec, and roll acceleration, rad/see^ , respec- 
tively 

amplitude of roll rate and yaw rate, respectively, rad/ sec 

amplitude of rocking-mode roll rate in total roll response 
to yaw -mode excitation (fig. 15(b)), rad/sec 

yaw-mode amplitude ratio of roll-to-yaw rate 
1 

\5 


XZ 


cos 6^^ sin 8^^ 
sp sp. 
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R 1 
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r,r 
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vertical reaction loads on main landing wheels and nose- 
wheel, respectively , in center-of-gravity determination 
(figs. 1(a) and 1(b)), N 

spring configuration in roll tests with springs normal to 
moment arm extending from spring attachment point on 
airplane to pivot of oscillation (fig. 6(c)) 

yaw rate, rad/sec, and yaw acceleration, rad/sec* , respec- 
tively 

time, see 


W 


weight , N 



sin 5 



ec 


gravity-force moment of experimental component , N-m 


W , weight of the aircraft with test-setup components removed , 

a/c 



weight of the aircraft with the weight of the test-setup 
components included , N 



weight of suspension beam (fig , 3) , N 


weight of ballast (fig. 3), N 
weight of outrigger (fig . 3) , N 


w 


X,Y,Z 



X ,Y ,Z 
s s s 



loading weight in the determination of the vertical center- 
of-gravity position using the single-point suspension 
method, N 

axes of the airplane through its center of gravity (fig. 17) 

longitudinal and vertical reference axes , respectively , 
through the suspension bolt (fig . 3) 

longitudinal principal axis of the airplane 


axes of the airplane (with test-setup components included) 
passing through the center of gravity located on the line 
of suspension (fig . 17) 

horizontal and vertical distances , respectively , of the 
center of gravity of the suspension beam relative to the 
X''Z-axes (fig. 3), m 
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X ,y ,z 
s s s 


X'' ,z^ 
w w 
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Ax, Az 
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St W 


, m 


I y I 

l^stl 


horizontal and vertical distances, respectively , of the 
center of gravity of the leveling ballast relative to the 
X^'Z^'-axes (fig . 3) , m 

horizontal and vertical distances, respectively, of the 
center of gravity of the spring attachment outrigger 
relative to the -axes (fig. 3) , m 

horizontal distance from line of suspension to the body 
reference station (fig . 3) , m 

horizontal, lateral, and vertical distances , respectively, 

of the center of gravity of an experimental component 

from the X Y Z -axes , m 
s s s 

horizontal and vertical distances, respectively, of the 
loading weight , w, relative to the X""z"-axes (figs. 1(c) 
and 3) , m 

distances parallel and normal, respectively, to the center 
of gravity of the aircraft from the datum used in ground- 
based determinations of the center of gravity (fig. 1(b)) , 
m 

shift in X and ~z, respectively, due to the removal of test- 
setup components , m 

lateral displacement , m 


sway-mode amplification factor 



z 


r 


z 


r 



amplitude of oscillatory lateral displacement due to sway 
mode, m 

vertical location of the center of gravity of the aircraft 
relative to the X-axis (fig. 1(c)), m 

vertical distance from the suspension bolt to the horizontal 
reference body axis defined by the reference rivets 
(fig. 3) , m 

vertical location of the center of gravity of the aircraft to 
the horizontal reference body axis defined by the refer- 
ence rivets , m 

value of aircraft with test-setup components included , m 
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inclination angle, relative to the horizontal, of the line from 
the point of attachment of the spring at the aircraft test set 
up to the pivot of oscillation (in pitch and roll tests) , rad 

change in 8^^ due to oscillatory perturbations , rad 

amplitude of A8j^, rad 

acceleration of A8j^, rad/sec^ 

inclination angle of spring plane (fig. 11) , rad 

inclination angle of inverted-pendulum arm of component 
weight, W (fig. 4), rad 

change in 8^ due to oscillatory perturbations , rad 

inclination of principal X-axis relative to aircraft X-axis 
(fig . 11) , rad 

pitch attitude , rad 

roll angle , rad 


roll angles of arms of double-pendulum concept of the 
single-point suspension (fig . 13) , rad 

amplitudes of sway-mode roll and yaw motions , respec- 
tively, rad 

amplitudes of rocking-mode roll and yaw motions , respec- 
tively, rad 


sway-mode amplitude ratios of roll and yaw motions of the 
aircraft , respectively , relative to lateral motion 


roll-to-yaw amplitude ratios of the aircraft for roll-mode 
and yaw-mode transient responses , respectively 

yaw angle, rad 
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to 


frequency of oscillation , rad/ sec 


to. 


frequency of oscillation of an uncoupled mode of oscillation , 
rad/sec 


CO ,C0 ,00, 

y <P t 




frequency of oscillation of sway-mode , roll-mode , and yaw- 
mode transient responses , respectively , rad/sec 

frequency of oscillation of simple, uncoupled sway-mode, 
roll-mode , and yaw-mode transient responses , respec- 
tively, rad/sec 


Aoo^ = 00^ - cOq , rad/sec 


Subscript: 

ec 


experimental component 


CENTER OF GRAVITY 


The horizontal position of the center of gravity of an airplane can be readily 
determined at the time the airplane is weighed . The vertical position can be deter- 
mined in several ways . The most commonly used techniques are discussed in this 
section and shown in figures 1(a) to 1(c) . 


Use of Ground-Based Loading Equipment To Determine Horizontal 
and Vertical Position of the Center of Gravity 

The weight of an airplane and the horizontal position of the center of gravity can 
be obtained by placing the airplane on loading cells or platform scales with the refer- 
ence X-axis horizontal and measuring the reaction forces on the two main landing 
wheels , Rj^ , and the nosewheel , R^^ , and the horizontal distance from the cen- 
ter of the main wheel bearings to the center of the nosewheel bearing, d (fig. 1(a)) , 
Thus 


w = Rn-Rm 


( 1 ) 


and the distance from the main wheel bearings to the center of gravity , x , is ob- 
tained from 


- 


( 2 ) 


The precision of the results depends on the minimum count resolution of the scales , 


9 



the accuracy of the measurement of d , and the precision with which the position of 
the main wheel bearings is known with respect to the body station . 

The vertical position of the center of gravity relative to a reference datum , as 
well as the horizontal position and the airplane weight , can be obtained by arranging 
the setup shown in figure 1(a) so that it allows the airplane to be tilted in the pitch 




(b) Weight and vertical and horizontal position of center of gravity . 



Reference X-axis 


(c) Vertical position of center of gravity . 

Figure 1. Experimental techniques for determining weight and 
center-of-gravity positions . 
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plane and then measuring the tilt angle (fig. 1(b)) . The reaction forces on the 
scales, and Rj^, and the horizontal distance, d, vary as the airplane is tilted. 

The tilt angle is obtained by means of an inclinometer or a transit. 

As shown in figure 1(b) , the center of the main wheel bearings is the datum 
point and a datum line is extended from the datum point normal to the reference 
X-axis of the airplane . Thus 

Rj^d = w(x cos Q sin g) (3) 


where x is the distance parallel to the reference X-axis from the datum line to the 
center of gravity , and z is the distance parallel to the datum line from the datum 
point to the center of gravity . Then from equation (3) , 


X - z tan 0 


V 

W cos 0 


(4) 


A graphical solution for x and z in equa- 
tion (4) , using several inclination angles , 
is obtained as shown in the adjacent sketch. 

Figures 2(a) and 2 (b) show the exper- 
imental setup , equations ,-and results of a 
test in which the weight and vertical and 
horizontal positions of the center of gravity 
were determined for the Fairey Delta 2 air- 
plane without fuel (ref. 10) . In this exper- 
iment , a point on the body was selected as 
the datum point. 




Figure 2. Example of the experimental determination of the horizontal and 
vertical position of the center of gravity of the Fairey Delta 2 airplane by 
tilting the airplane on weighing scales. Data are from reference 10. 
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(b) Horizontal and vertical position of center of gravity. 
Figure 2. Concluded. 


When using this technique , precautions must be taken to insure the rigidity of 
the support points relative to the body and the precision of the measurement of the 
distances and attitude angles . Otherwise , the results may be unreliable . If flexure 
at the loading points cannot be eliminated , other means are necessary to determine 
the center of gravity . 


Use of Single-Point Suspension To Determine the 
Vertical Position of the Center of Gravity 

The vertical position of the center of gravity , as well as a refinement of the 
horizontal position , can be determined as byproducts of the use of a single-point sus- 
pension double-pendulum rig to obtain the inclination of the principal axis and the 
moment of inertia about the Z-axis. In this technique, the reference X- and Y-axes 
of the suspended airplane (fig . 1(c) ) are leveled by the addition of ballast (normally 
lead shot) . While the airplane is being leveled , which is usually done with a clin- 
ometer , care must be taken not to cause pitch attitude changes due to the load in- 
duced by the clinometer operator . After the airplane is level , a weight box of known 
weight is suspended from the front of the airplane with the point of suspension known 
accurately both horizontally and vertically with regard to the pivot of suspension. 
Known weights are added to the weight box, and the stabilized vehicle’s attitude for 
each weight loading is measured . For each weight loading (which includes the 
weight box) and corresponding pitch attitude , the vertical position of the center of 
gravity relative to the pivot, without the weight loading, is determined as follows: 


w(x'' cos 0 - z'' sin 0^ 
— \ w w / 

W sin 9 


W \^tan 0 ^w^ 


(5) 


where w is the weight loading , and W is the suspended aircraft weight with ex- 
perimental gear but without loading . 
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To insure consistent results with this technique , the following precautions 
should be observed: 

(1) The pivot point should be as free from friction as possible. 

(2) The vertical distance from the pivot to the suspension point of the weight 
box should be as short as practical in relation to the horizontal distance to insure 
the maximum tilt angle for each weight condition . 

(3) The sling lines of the suspension should converge and be fastened to a col- 
lecting fixture , which itself contains the eye for attachment to the suspension hook . 

(4) The measurement of the tilt angle for each loading should be as accurate as 
possible. Good results have been obtained by using a transit to sight the vehicle’s 
attitude by reading two rigid tape measures attached at reference points along and 
normal to the reference X-axis . The tapes are read at each stabilized tilt position , 
and the difference in the readings divided by the distance between them (measured 
along the X-axis) provides the tangent of the tilt angle . Using an inclinometer to 
measure the tilt angle has not been satisfactory . 

(5) To minimize the effects of hysteresis , the loadings should be removed in 
reverse order of increasing load application , and the average reading for each tape 
at each load condition should be obtained . 

Figure 3(a) shows the suspended setup used at the NASA Flight Research Center 
with the HL-10 lifting body vehicle. Figure 3(b) shows the vertical position of 


Wa,c = 23, 343.75 N 
W|j = 733.95 N 
Wb3| = 142.34 N 
W^j./outrigger = 44.48 N | 
(Wg/c) = 24,309.00 N 

x^ef = -0.359 tn _[ 

l^ = 2.964 m z' = 0.884 m 

x' = -0.419 m z' = 0.162 m 
b b 

xLi = -0.089 m zL, = 1.892 m 
bal bal 

x'|. = -0.051m z'|. = 1.892 m 

x'=3.110m z'= 1.337 m 

w w 


l/-Suspension line 



-Front tape for measuring 
angular deflection with 
transit 

Center of gravity of outrigger 


Rear tape for measuring 
Reference angular deflection with 
station transit 
3.632 m 

Center of gravity of ballast 


(a) Test setup . 

Figure 3. Experimental setup and results obtained for the vertical position of 
the center of gravity of the HL-10 lifting body vehicle by using the single-point 
suspension method. 
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(b) Vertical position of the center of gravity with experimental gear included. 


Figure 3. Concluded. 


the center of gravity of the vehicle when weighted with the leveling ballast and the 
experimental gear as determined for each tilt loading , as well as the average value 
from all loadings . The average value is within ±0.005 meter of the spread of the test 
results. Appendix A summarizes the data, calculations, and results of determining 
the vertical and horizontal positions of the center of gravity of the vehicle. 


PITCHING AND ROLLING MOMENTS OF INERTIA 


The determination of pitching and rolling moments of inertia with the airplane 
pivoted on knife edges involves perturbations of the gravity-force moments (weight 
moments) , restoring spring moments , and corrections for the apparent additional 
mass effect of the ambient air mass set in motion by the motions of the airplane . For- 
mulas for the correction of apparent additional mass effects on the test results are 
given in reference 11 . 

If the nonlinear effects of perturbed gravity-force moments and spring moments 
have been minimized to the point of insignificance , the equations of motion for a 
knife-edge setup of a compound pendulum for the determination of pitching-moment 
inertia are as follows: 



+ 


w , 

a/c „2 
g a/c 



sin 5, 


W> 


A6, 


ec 


( 6 ) 


Letting 


AS^ = sin to^t 


where co^ is the natural frequency of the oscillatory response for the setup , and 
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transposing equation (6) results in the expression 


K 


W 2 ta 7 sin 6 
W 


W, 


ec 


W 


00 


a/c 

S 




( 7 ) 


The equation for rolling moment of inertia has the same format as equation (7) . 


The pitching moment of inertia about the pivot axis of the experimental rig 


components , ^ 


0 


g W 


) ’ 


, and, in the case of rolling-moment inertia, 


, are either calculated or, if feasible, determined experimentally . 


The apparent additional mass is calculated by using the formulas in reference 11. 


It remains to determine the effect of the amplitude of oscillatory perturbations , 

|A6j^|, in conjunction with the position of the airplane’s center of gravity relative to 

the pivot and also the spring configurations , on the degree of nonlinearity in the 
responses to the perturbations . 


Gravity-Force Moments 

In a rolling-moment-of-inertia setup , the weight vector for the unperturbed 
condition normally lies in the vertical plane that extends through the pivot . During 
roll oscillations , the change in weight moments is accounted for by the follov^ing 
equation, which is for a simple inverted pendulum (see sketch) and shows the 
change in weight moments to be a linear function of roll displacement: 


Y 



/ 


/ 


W 





z 


AM^ = sin A8^ 




( 8 ) 


In a pitching-moment-of-inertia setup , the weight vector for the unperturbed 
condition is forward or aft of the pivot and the equation takes on the following 
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nonlinear format: 



(9a) 


(9b) 


In equation (9b) , is the ratio of the nonlinear to the linear change in weight 
moment for the condition when the inverted pendulum is inclined at an angle , 8^ , 
under steady-state conditions . The ratio is a function of both the magnitude 
and the sign of the perturbation angle , A8^, as shown in figure 4. As the figure 
shows , the effects of nonlinearity can be significant with increasing in one 

direction. However, during oscillatory motion, if A6^ has the same positive and 
negative amplitudes , the algebraic subtraction of equation (9a) results in the 



Figure 4. Ratio of nonlinear to linear increment of change in weight moment 
as a function of steady -state and perturbed inclinations 8^^ and A8^, re- 
spectively (eq. (9b)) , 
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double- amplitude expression 2AM^ = 2W2^ sin 6^ | AS^ j . Thus , the nonlinearity 

would cancel out. On the other hand, the algebraic addition of equation (9a) results 

in a net = 2W2^ sin 8 ^|e^| instead of zero, indicating that there is a tendency 

for the zero datum of the time history of the pitching-moment perturbations to shift . 

It should be noted that when approaches 0° , equation (9a) reduces to the non- 
linear expression 

8w=0 


Spring Moments 

A number of options are available in the layout of the springs in the pitching- and 
rolling-moment-of-inertia test setups , For pitch setups , the single spring (or its 
equivalent) in a simple spring installation may be forward or rearward of the pitch 
pivot . For roll setups , identical sets of springs are usually attached to the vehicle 
on each side of the pivot below the spring mounting points of the experimental setup 
to provide a stabilizing restoring couple during roll perturbations . Two sets of 
springs can also be used in pitch setups . 

Regardless of the spring setup, the springs must be carefully calibrated and 
pretensioned to make sure that they operate within their linear range during pertur- 
bations . A typical calibration curve is shown in figure 5 . 



Figure 5. Typical calibration curve of a tension spring 
with looped ends, k = 7544 N/m. 
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Springs normal to the moment arm at point of spring attachment. — For single 
spring configurations , used in pitching-moment-of-inertia investigations , with the 
springs normal to the moment arm that extends from the spring attachment point to 
the pivot (figs. 6(a) and 6(b)) , the perturbated spring moment can be represented 
by the following equation: 


AMg = -kg2A8^(l + Eg) (11) 

where k is the linear spring constant, £. is the length of the moment arm, and E 

X S 

is the ratio of the nonlinear increment in the change in spring moment to the linear 
change. For the configuration identified as P^l in figure 6(a) , 



(12a) 


For the configuration identified as P^2 in figure 6(b) , 



(12b) 


Equations (12a) and (12b) and the plotted ratios in figures 6(a) and 6(b) show that 
the only difference in configurations P^l and P^2 is the change in the sign of 

A8^. For both configurations , the direction of the perturbation, A8j^, affects the 

magnitude , but not the sign , of E . 

s 

For amplitudes of A8^ of 2° or less , E is negligible (-0.001 or less) . Even at 

JL S 

A8^ = 4°, E is only approximately -0.005. However, for large amplitudes of A8. , 

X S X 

10° for example, the ratio is significant. 

Figure 6 (c) shows spring configurations for pitching- and rolling-moment-of- 
inertia investigations in which the configuration has identical springs on each side of 
the pivot and identical moment arms normal to the springs from the spring attachment 
points to the pivot to provide a restoring couple . In these configurations , the change 
in spring moment due to perturbations can be determined from the expression 


AM3 = -2M^jA6j(1 + eJ 


(13) 
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(a) Spring configuration P^l for pitch. 


Figure 6. Ratio of nonlinear to linear increment of change in spring moment 

h 

as a function of A6 and ^ . Spring normal to moment arm at point of 


attachment to airplane. 
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where E is obtained from figure 6 (c) or the following equation: 
s 


E = 


s 




( 14 ) 


It is apparent that E is a function of both 
sign , of A8 ^ . 



and the magnitude , but not the 


Although configurations with springs normal to the moment arm at the spring 
attachment point are generally considered to provide linear spring moment charac- 
teristics over a large range of perturbations , figure 6 indicates that this is not nec- 
essarily true. As shown below , the nonlinear effects are negligible over a larger 
perturbation range for some vertical-spring configurations . 


Vertical springs . — For single vertical-spring configurations used in pitching- 
moment-of-inertia studies , the change in the spring moment due to perturbations is 
determined from the following expression in which it is desired to minimize the non- 
linear error ratio, E": 

s 


AM = -kC? A6, (l + E ( 15 ) 

s Ij^ 1 \ s/ 

Two configurations , P^l and P^2 , of single vertical springs mounted forward of 

the pivot (and their mathematical equivalents with the springs rearward of the pivot) 
are illustrated in the following sketches , which also indicate the positive direction of 
the quantities shown: 




P^l configuration 




P^2 configuration 
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For configuration P 1, the expression for E is 

V s 


^2 ■? 


\2 


) \2 

v) 






tan 

(A6^) tan^ 6^ 


E =- 
s 




(16) 


Since the P^2 configuration is an inversion of the P^l configuration , changing the 
signs of 8. and A8^ in the equation for P 1 results in the following equation for 

•n rt . t t V 


E 


s 


r 


2 3 ,.j. 




'*^2V 




\V ’ 


tan 8j^ ^ 
(A8^) tan^ 6^ 



(17) 


It should be noted that when 5^^ = 0, the equations reduce to equations (12a) and 

(12b) , which are for configurations in which the springs are normal to the moment 
arm at the spring attachment point . 

The nonlinear increment ratio, E , is plotted in figures 7(a) to 7(h) as a ■ 

^2 

function of tan 8^^ and for a series of values of ^ — for configuration P^l- 

^h 

The values of E for configuration P 2 can be obtained from figures 7(a) to 7(h) 

S V 

by reversing the signs of tan 8. and A8^. The values of E in the figure for 

X X s 
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tan = 0 are the values shown in figures 6(a) and 6(b) for the configuration in 
which the springs are normal to the moment arm that extends from the spring attach- 

*2 

ment point to the pivot . Figure 7 shows that for each value of -s — there is a trend 

ih 

toward a crossover point , which is also a function of tan 8. , where E = 0. These 

crossover points represent the optimum design conditions for the single vertical- 
spring configuration, and when so configured, the single vertical- spring setup is 
superior to the spring configurations shown in figures 6(a) and 6 (b) insofar as the 
effects of nonlinearity are concerned . 



(a) y^=0.2. 


Figure 7. Ratio of nonlinear to linear increment of change 
in single vertical-spring configuration P 1 as a function 

1^1 
of tan 8 y A8^, and j — for several values of y— . 


24 



■ 1.6 - 1.4 - 1.2 - 1.0 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 

tan 6j 


(C) 



0 . 8 . 


Figure 7. Continued. 



















(h) 



4.0. 


Figure 7. Concluded. 


Figure 7 also shows that if E is to be held within a prescribed limit (for ex- 

s 

ample, less than 0.005) , the permissible magnitude of the perturbation is governed 
by 6. and -p— . 

ih 


When the vertical springs are mounted symmetrically about a vertical plane of 

symmetry that passes through the pivot axis , as in rolling-moment-of-inertia setups , 

E in the spring moment equation 
s 


AM =-2k«2 A5,('1 + eJ 
s 1, 1\ s/ 


is defined by the expression 



(18) 


( 19 ) 
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where , as shown in figure 8, 



( 20 ) 


It is apparent from equation (19) that E can be eliminated by making B equal to 

s 

zero . This can be accomplished by making the length ©Qual to twice the length 
of gg (eq. (20)); that is , by arranging 9 .^ so that it is bisected by the horizontal 
projection from the pivot. 


The nomograph in figure 8 can be used to obtain values of E for different 

o 

values of A8. for any two-spring couple configuration. For a limited prescribed 

^ g 

^0 

value of E and a fixed value of s— , the nomograph can be used to obtain the 

S M .4 

h go 

limiting values of A6. for different values of . For a prescribed limiting value 

of E at a set value of «— > it can be used to obtain the value of p— for a selec- 

- ^h h 

ted upper limit of AS^^ . For example , for the configuration shown in the upper left 

of the figure , in which the springs are below the attachment points , let E ^ 0 . 005 , 
g “ 

g — = 0.6, and A8j^ = ±4.5°. From the right side of the nomograph, find the inter- 

^h g2 

section of A6-, = ±4.5 and «— = 0.6 and draw a horizontal line to the left. From 

ih 

the left side of the nomograph, draw a vertical line from E = -0.005 and note the 

s 

value of B at the point where the lines intersect. In this example, B = ±0.8. Now 
if 



then 


from which 


± 0.8 = 0.6 - 




= -0.1 or 0.7 
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Positive 



0 .1 .2 .3 .4 .5 .6 




Figure 8. Variation of nonlinear increment ratio of two vertical 
springs, configured to provide a restoring couple, as a function 


of A8^, 



and 
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or , as shown in the following sketches , 


or 


Positive 




If the stipulated design conditions are applied to the configuration shown in the 
upper right of figure 8, in which the springs are above the attachment point , the 
results would be a mirror image of the preceding results . 


Practical Considerations in Setting Up the Experimental Installation 

Unless care is taken in every detail of an experimental setup for pitch and roll 
investigations , difficulties may be encountered . During one experiment to determine 
the pitching moment of inertia in which the aircraft was supported at the wing jack 
points and oscillated with the spring at the nose, the wing, which had been consid- 
ered rigid, flexed as the aircraft oscillated . This flexing caused the axis of rotation 
to shift forward and downward from the line through the jack points. 

A common error is using flexible cables to tie back the springs . Flexible cables 
should not be used, since they can constitute springs in series with the actual 
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springs , resulting in a softer spring system than intended or accounted for. Using 
rods for tiebacks alleviates this problem. It should also be noted that on some air- 
craft the rear portion of the fuselage is relatively flexible , so that attaching the 
spring to it alters the effective spring constant . The same problem can arise if the 
springs are attached to the wings . Thus , the springs should be attached to the 
cradle, or, if they are attached to parts of the structure that are suspected of being 
flexible, bracing should be set up from the cradle to the spring attachment point. 

Inertia characteristics are usually determined for no-fuel conditions because the 
sloshing of the fuel tends to introduce a beat action into the transient oscillatory re- 
sponse after excitation. When inertia determination is attempted with fuel on board, 
the difference between the oscillatory mode of the sloshing fuel and that of the experi- 
mental setup should be as large as practical to make it possible to isolate the beat 
action . 

Since the experimental setups for determining the pitching and rolling moments 
of inertia can involve substantial effects of such experimental accessories as cradles 
(fig . 9) , the vertical and horizontal positions of the center of gravity for the clean 



liliil; 










vehicle and the experimental components must be known . The center-of-gravity 
positions for the clean vehicle can easily be determined experimentally by following 
the procedures discussed earlier . The horizontal and vertical positions of the center 
of gravity of the experimental accessories can be calculated and verified , at least in 
part, by experiment. 


DETERMINING THE PRODUCT OF INERTIA AND YAWING MOMENT OF INERTIA 
BY THE SINGLE -POINT SUSPENSION DOUBLE-PENDULUM METHOD 


Until the single-point suspension double-pendulum method described in refer- 
ence 7 was introduced, the experimental determination of the product of inertia was 
difficult and usually unsatisfactory , In the single-point suspension double-pendulum 
method, which provides more accurate results, the airplane is suspended by a sling 
from a point at the end of a suspension line (figs. 10(a) and 10(b)) . The suspension 


(a) Springs oriented crosswise. 


Figure 10. Experimental setup for determining the yawing moment of 
inertia, the product of inertia, and the inclination of the principal axis 
using the single-point suspension method. 







(b) Springs oriented lengthwise . E- 14681 


Figure 10. Concluded. 


line is , in effect , a torsion-free cable . The other end of the line is attached to a 
crane hook. The sling system itself is a rigid assembly. Yaw and roll restraint is 
provided by two sets of springs whose lines of action lie in a common plane 
(fig. 11(a)) . The springs are selected and arranged to provide a spring couple 
vector in the XZ -plane of symmetry, which acts normal to the spring plane. One 
component of the vector provides restraint in yaw; the other provides restraint in 
roll (fig. 11(b)). 

The crosswise spring configuration shown in figure 10(a) was used in the study 
of reference 7 and was given extensive theoretical treatment in reference 8 . The 
longitudinal spring configuration shown in figures 10(b) and 11(a) , which is the 
arrangement used at the NASA Flight Research Center , is considered in this report . 
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(a) Test setup. 



(b) Resolution of spring couple vector. 


Figure 11. Test setup and resolution of spring couple vector in single- 
point suspension method for determining the yawing moment of inertia, 
the product of inertia, and the inclination of the principal axis . 
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Objective and Principle of the Method 

The objective of the single-point suspension method is to determine the inclina- 
tion of the spring plane at which the roll-to-yaw ratio of the yaw mode of oscillation 
is zero . After this inclination and the corresponding frequency of the yaw mode of 
oscillation have been determined, the yawing moment of inertia, the inclination of 
the principal axis, and the product of inertia can be determined. 

In principle , as outlined in reference 7 , the method considers a pure yaw-mode 

response to a yaw excitation. With the spring plane at some one inclination, 8_„ 

sp 

(fig. 11(b)) , a yaw-mode excitation of the aircraft results in a certain amount of 

oscillatory roll response as well as yaw response . The restoring spring couple 

vector, C (fig. 11(b)) , produces a restraining moment on the yaw motion and a 
s 

rolling moment that either reinforces or restrains the rolling moment, depending on 
the inclination of the vector . The yaw-mode transient responses in roll and yaw , 
following yaw excitation, are represented by the following equations: 


‘xP ■ 'xz' = 


I„r - = -C cos 6 

Z XZ^ s sp 


(21a) 

(21b) 


At some one value of 8 

sp 

condition cp = 0 and thus 


the rolling motion of the yaw mode will be zero . 
p = 0 , and the equations reduce to 


At this 


->xzf = 


I„f = ” C cos ^6 
Z s I spj^^p 


(22a) 

(22b) 


Hence , when cp = 0 and thus = 0 , 

r\^ 


U 


tan (8 ) -j 

\ sp/„_„ I 


XZ 

cp=0 h 


(23) 


and 

^XZ " h 

Since the inclination of the principal pitch axis for a symmetrical airplane is given 
by the expression 

tan 2e = , T (25) 

Z X 
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where £ is the inclination of the principal axis relative to the horizontal reference 
X-axis , substituting equation (24) into equation (25) gives 

2I„ tan (8 \ 

Z V sp ) 

tan 2e = ^ ^ ^ (26) 


The yawing moment of inertia can readily be determined as a byproduct of 


this method . 

ril 


Using equation (22b) and considering 5 to be the angle at which 


sp 




0 , 


“C cos 
s 




(®sp) 


(p=0 


(27) 


With the springs configured longitudinally, as in figure 11(a) , 
is represented by the expression 


for 


l_9i 


= 0 


c = 

S 



COS 



(28) 


and f is represented by -vj/ / co . V . Substituting these expressions into equa- 
tion (27) results in ^ 


4k(8')' cos 


(29) 


Response Characteristics of the Actual System 

A simple yaw and roll motion response to a yaw-mode excitation at the yaw-mode 
frequency would be ideal for the single-point suspension method . Actually , the sys- 
tem allows the aircraft to rotate about all three axes and to translate along the pitch 
and roll axes . The proper application of excitation can eliminate the longitudinal 
motion in a lateral excitation; however, the yaw mode, the rocking mode, and the 
sway mode of oscillation (fig. 12) must be accounted for. 

The yaw mode is primarily a yawing motion at the yaw -mode frequency, with 
some coupled roll motion and a small amount of sidewise motion . The yaw mode is 
required for the determination of 1^^ ^2 ‘ 

The rocking mode , ideally , is a rolling motion at rocking-mode frequency with- 
out sidewise or yaw motion , Actually , it may contain a small amount of sidewise and 
yaw motion . The unintentional excitation of this mode during the excitation of the 
yaw mode produces a significant and unwanted interference with the roll response of 
the yaw mode in recorded time histories . This interference can be accounted for in 
the analysis of the response data. 


37 




Yaw mode Rocking mode Sway mode 


Figure 12. Rolling motions of the three natural modes of oscillation to be 
accounted for in responses to yaw-mode excitation in the single-point 
suspension method for determining the product of inertia and the inclina- 
tion of the principal axis. Motions are exaggerated for clarity . 


The sway mode iji primarily a sidewise motion at the sway-mode frequency , 
where it usually contains some roll motion and a negligible amount of yaw motion . 
The interference of the sway-mode roll motion with the yaw -mode roll motion can be 
made negligible by proper design of the suspension rig . 

Since the roll-to-yaw ratio and 
the frequency of oscillation of the yaw 
mode are the desired response char- 
acteristics , and should be accounted 
for quantitatively as a function of the 
inclination of the spring plane , it is 
essential to obtain a mathematical in- 
sight into the interacting effects of 
the three modes of oscillation on the 
frequency of oscillation and amplitude 
ratio of the yaw mode . This insight 
enhances the minimization and isola- 
tion of the interference of the rocking 
and sway modes with the yaw -mode 
characteristics . 

Equations of motion. — The follow- 
ing spring forces and moments , which 
are used in the equations of motion, 
are based on a small-perturbation 
analysis of the longitudinal configura- 
tion of the springs and the schematics Figure 13. Schematic diagram of small- 
of the suspension system shown in perturbation motions and gravitation force 

figures 11(a) and 13. Experience has in single-point suspension system. 
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shown that since the distance from the spring attachment point at the airplane to the 
spring stanchion (distance a in fig. 11(a)) is of the order of 3,5 meters to 
4.5 meters, the spring side forces are normally negligible in comparison with the 
side component of the gravitational force . Similarly , the spring moments due to the 
spring side forces are negligible in comparison with the spring moments due to 
angular motions . Thus , for the purposes of the present study, 

Y « 0 (30) 

s 


Ls ^ ®sp ^sp)'*' ■ (^t 

Ng ^ Sgp sin 6^p)9i - ( cos^ 6^^) (32) 


where 


k, = 4k(l!9‘ 


(33) 


The contributions of the gravitational and inertial forces and moments to the equa- 
tions of motion are obtained by using Lagrangian equations and figure 13 . The 
following equations of motion can then be obtained by combining the spring with the 
gravitational and inertial contributions: 


I 


0 


/d^ + 


‘pi 


Y 



g 

[ 



s 


/ 



H'l 







-I„ D* 


w— 





-I 

XZ^ 


IzD^ 


0 


y 


Ns 






_ 






(34a) 

(34b) 

(34c) 


By substituting equations (30), (31), and (32) for Y^, and N^, respectively, 
and consolidating. 



(35 a) 
(35b) 
(35c) 
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where 


A = W 



= T-^ COS 8 sin 8 

^ ixz ®P 



2 


0 £ 


0 







2 

0 



8 


sp 


(36) 

(37) 

(38) 

(39) 

(40) 


Yaw-mode frequency . — From the matrix of equations (35) , the characteristic 
equation can readily be determined to be 

[“■ * (%);] [”■ * {“*);] [“■ ‘ (“y)J - [“■ * [■>■ • (“y)’„] - - (“t)*.] = » <«> 

A direct solution of equation (41) , for yaw-mode frequency , is complex. A simpler 
solution , described in reference 8 , considers a mode frequency on the basis of 


CC^ = + Ao)^ 


(42) 


Thus, to obtain the yaw-mode frequency, co^, from the characteristic equation , 
let - A . The resulting solution gives the following 

expression for A • 


XZ 


[(“♦)’ ■ p!]' 




X Z 




2 _ 


(“<.)’ - (%) 


(43) 


1 - 


g X 


■ (“<p)j] [(“t)' - ("y);’ 
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At the expense of a small error in A/oo . ^ , (co.^ can be substituted for 
the equation. V ¥/ \ ¥/q \ ¥/ 


in 


Equation (43) indicates that the yaw-mode frequency is affected to some extent 

by the uncoupled roll-mode frequency, , and the uncoupled sway-mode fre- 

\ 9/q 

quency , (oo V , except at the one condition pertinent to the experimental Investi- 

\ Y/q 

gation . When p^ equals (w . V , which is the condition required for the yaw -mode 
^ |cpl I 

amplitude ratio, ^ -' j- , to be zero (as shown later) , A/co . equals zero and the 

I Y 1^ V Y/ 

yaw-mode frequency is simply > which is the frequency of a simple uncoupled 

yaw motion. This is fortunate, since this uncoupled frequency, obtained experi- 
mentally, is required for the accurate determination of the yawing moment of inertia, 
which is essential in obtaining the product of inertia . 

Amplitude ratios . — The roll-to-yaw amplitude ratio of the yaw mode is obtained 
from equation (35b) after substituting equation (35a) for y in equation (35b) and 
letting equal . Thus , 


XZ 


[(“♦)’ - pl] 




l^ll 

l'^l¥ 


(“♦)' - (“<p) 


(44) 


1 - 




Equation (44) indicates that there should be a significant difference between the 
magnitude of /co , V and that of • In addition, the equation indicates that 

^ ^^0 icp.i 

when p^ equals ( . the amplitude ratio j | is zero . At this condition , with 

the corresponding inclination angle, 6 , of the spring plane and the uncoupled 

sp 

yaw-mode frequency, , known from experimental data, the yawing moment of 

inertia is obtained from equation (40) and the product of inertia is obtained by 
equating equations (37) and (40) . Thus, 


cos^ 




(®sp) 


cp=0 


(45) 


(“-I'): 
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and 





( 24 ) 


The roll-to-yaw amplitude ratio of the rocking mode is obtained from equa- 
tion (35c) with equal to , Thus, 


,1 [(“.); - (%)’; 
rn<p 'xz [pi-(<»,p)’] 


(46) 


Equation (46) indicates that the rolling motion of the rocking mode can be elimin- 
ated by making equal to ^ to avoid possible interference of this rolling 

motion with the roll response of the yaw mode . However , this is not feasible , since 
it conflicts with the yaw-mode requirement of equation (44) for a good amount of 
separation between (“cp)^ ' ^or the critical condition of the yaw mode 

= 0 , at which condition p^ equals /co . , equation (46) reduces to 

I '1^ l\|r _ 


9 




9 


XZ 


(47) 


At this critical condition , the rocking mode has a high roll-to-yaw ratio that is deter- 
mined by the inertial characteristics alone . The amplitude of the rocking-mode roll 

motion » 1 9 I > is dependent upon the amplitude of the rocking-mode yaw motion , 

11^ 

j \1/ , unintentionally induced during the excitation of the yaw mode . 


For the sway mode , the amplitude ratios 


9 


and 




are the ratios of 


rfly Fy|y 

interest in obtaining insight into sway-mode interference with the yaw-mode ampli- 


tude ratio. 


9 


1^' 


The following expression for 




-y 


can be obtained from equa- 


tions (35) with approximately equal to 


9i 


"X 




( 48 ) 


Substituting equations (36), (38) , and (39) for A, > and , 
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respectively , results in 


'I’ll 


ro^'O 


k. sin^ 6 gl„ 
t sp _ °X 


W 




wc' 


(49) 


Since the amplitude of the sway-mode rolling motion, |<Pj^ I . is primarily dependent 

y 

upon both Kq + £ and the amplitude of the sway-mode sidewise motion, pY ly > the 

minimization of this unwanted motion is dependent upon the care taken to avoid 
exciting the sway mode during the excitation of the yaw mode as well as upon the 
magnitude of 2^ + , 

The sway-mode amplitude ratio , 1^77 , can also be obtained from equa- 


|-y| 


y 


tions (35b) and (35c) with approximately equal to -/to V : 

\ y/Q 


JlA - 

-y\y' 


A^xz\ 


Pi - ( 

CO 

yloj 


\ ^xh) 

>t); - 



(50) 


|9i| 

For the critical condition of the yaw mode , ^ ^ | = 0, at which condition p^ = 


U I 

|-y|. 


(AI 


XZ 




(51) 


or 


hH . 

-y|y' 


I 


XZ 




k. sin^ 6 gl 
t ^ ® 


X 


C' 

1 

w 4 

*^0 




(52a) 
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(52b) 


I « I , 'xz 

|-y|~ 'z 


1 



For this condition, the sway-mode yaw amplitude is negligible. It can also be con- 
sidered to be negligible for conditions other than the critical condition . 

As indicated by the preceding discussion , the sway-mode yawing motion is neg- 
ligible and will not interfere with obtaining yaw-mode roll-to-yaw ratios from experi- 
mental time histories . The sway-mode rolling motion can also be reduced to negligi- 
bility by proper suspension system design . The minimization of the rocking-mode 
rolling and yawing motions depends primarily upon the technique used to excite the 
yaw-mode motions . When rocking-mode rolling and yawing motions are present in 
experimental time histories , it is relatively simple to isolate them to obtain the 
desired yaw-mode data. 


Minimization of the sway mode. — To minimize the effects of the sway-mode roll- 
ing motion on the total rolling motion in attempting to obtain yaw-mode responses , 


the sway-mode amplitude ratio, 
be minimized, since 



and the amplification factor , 



must 


fl\ _ l^ll I y I 

I ^st|y I ^ly l^stl 


(53) 


Reference 8 

prevent the 
ing motion . 


suggests that 



should be less than 0 . 164 radian per meter to 


sway-mode rolling motion from having signiflcant effect on the total roll- 
Using this as a criterion and applying it to equation (48) , 


0.164 > 




(54) 


By substituting equations (36) , (38) , and (39) for A, ( “ > and ( to , respec 

V y/Q V T/q 

tively, equation (54) can be used to solve for + 8 ^^ for a selected value of 8 ^^; 


8j +8^> 6.1 - 


k. sin^ 8 ^ 

_t S£ , g W 

z' 8 :: 

z' 


(55) 


The first term in equation (55) is dominant . At the NASA Flight Research Center , 
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inertia experiments are normally conducted with a value of of approximately 

7 . 5 meters . This is in accord with the criterion . 

The expression for the amplification factor j - ^^ 4 is obtained by introducing a 

r stl 

sinusoidal exciting side force , Yp , applied at the yaw -mode frequency , into equa- 
tion (35a) and considering the yaw angle to be equal to zero in the rolling-moment 
equation (eq. (35b) ) . The exciting force is considered to be so positioned as not to 
excite the rocking mode . Thus equation (35a) becomes 

Acp ^ ^ mD^ + y = Yp (56 a) 

and equation (35b) becomes 




= 0 


(56b) 


With equal to , where is the exciting frequency, and 

the following equation can be obtained: 


F 

Ygt fo IT ’ 
y 



(57a) 


(57b) 


It is apparent that the attenuation of j-^ - j requires ( to , V to be much greater 

than • This condition is readily satisfied, since is normally greater 

than 12 and f to V , which is equal to -S- , can be made desirably smaller by making 
V y/Q 9 ' 
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Kq as long as possible or 2 ^ as short as possible . 

The effects of the sway mode on total rolling motion can be reduced further by 
minimizing the excitation of y^^ . This can be done by gradually building up a 

sinusoidal yaw motion at the yaw-mode frequency with an exciting force that is kept 
small by placing it as far forward or rearward of the center of gravity as possible . 

Minimization and isolation of the rocking-mode rolling motions. — If the vertical 
positioning of the exciting force is improper during the excitation of the yaw mode , 
the exciting force may also excite the rocking mode . The factors that affect the 
response of the rocking-mode rolling motion to a lateral force that excites the yaw 
mode at the yaw-mode natural frequency can be obtained from the following approxi- 
mation of the rolling-moment equation: 

where 


L„ = k CD , 
F 9^ St 


(59) 


With equal to , where co^ is the exciting frequency, the following ex- 

pression can be obtained from equation (58): 



which reduces to 




(60) 


Equation (60) indicates that the rolling motion of the rocking mode can be atten- 

(co )2 

uated by minimizing 9 , and also by making — -x— greater than 2. Since 

((o V 
V Wn 
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cp is equal to tj- — , the minimization of cp depends upon keeping the yaw-mode 
cp 

exciting force as small as possible and positioning this force properly . Attenuating 


the rolling motion of the rocking mode by making 




much greater than 2 may 


0 


not be desirable because of the possibility of exciting structural modes of vibration , 
which should be avoided . In addition, should be dissimilar to , 

since equation (44) indicates that a well-defined difference is necessary if the yaw 

l^ll ^ Ip I 


mode 


i I” I4, 


, or its equivalent , 


|r| 




, is to be determined . At the NASA Flight 


Research Center , it has been the practice to make 


(%) 


approximately equal to 


0 . 5 to avoid the excitation of structural modes of vibration as well as to provide a 
reasonable separation between the two frequencies . 


From the preceding discussions of the single-point suspension, it is obvious that 
the frequencies of the sway mode, rocking mode, and yaw mode should be separated 
by a significant amount by the design of the suspension system . Also , the excitation 
of the yaw mode shouldTje sinusoidal and performed in a way that does not excite the 
rocking mode. In addition, the damping of the responses should be negligible. 
Finally, if the suspension system is properly installed, the roll- rate and yaw-rate 
responses (the responses normally recorded) are in phase or 180° out of phase . 

When rocking-mode rolling and yawing motions cannot be avoided, and are there- 
fore present in the total response to yaw excitation , they can be isolated from the total 
motion in the time histories (provided that the sway-mode motions are insignificant) , 
and the yaw-mode roll-to-yaw ratio can be obtained by plotting the peak values of 
the roll rate against the peak values of the yaw rate . In essence , this constitutes 
plotting the bounding outline of a lissajous figure . (Lissajous figures are discussed 
in detail in reference 8 . ) As shown below , this method of identifying and eliminating 
the rocking-mode motion from time histories to obtain the yaw mode can be applied 
even if significant damping is present . 

Figure 14(a) shows time histories of roll-rate and yaw-rate responses to yaw- 
mode excitation in which the rocking-mode interference with the yaw-mode rolling 
motion was absent but damping was present . A cross plot of the peaks of roll rate 

. plotted against the peaks of yaw rate (to obtain the of the yaw mode) is shown 

I ^Ix); 

in figure 14(b) . Because of the absence of rocking- mode interference with the yaw- 
mode responses , the cross plot shows only a single line. Also, because of the damp- 
ing, there are several points along the line; if damping had not been present, the 
line would have been defined by only two points , one at each end . 
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Time 


p, deg/sec 



r, deg/sec 



(a) Time histories of transient responses. 



(b) Lissajous cross plot; |-2| = 0. 485. 

Figure 14 . Time histories and lissajous cross plot of a lightly 
damped transient response free of rocking- and sway-mode 
interference after yaw-mode excitation. 


When the rolling motion of the rocking mode interferes with the rolling motion of 
the yaw mode and damping is negligible , the time histories of roll- and yaw-rate 
responses are like those in figure 15(a) (from ref. 8) . A cross plot of the peaks of 
these roll-rate and yaw-rate responses results in a lateral spread of the points , as 
shown in figure 15(b) . Drawing a boundary outline of the points, as shown in this 
figure , produces a parallelogram that is the bounding outline of a lissajous figure . 
The mean line of the left and right sides of the parallelogram represents the response 
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of the yaw mode . The slope of this line is the amplitude ratio of the yaw mode 



for the tilt angle of the spring plane , 6 . The lateral distance from this mean line 

sp 

to the side of the parallelogram is the amplitude due to the rocking mode, |p|jp • 


p, deg/see 


r, deg/see 





(a) Time histories of transient responses (from ref. 8). 



(b) 


Lissajous cross plot; 


r\ 


'I' 


0.078. 


Figure 15. Time histories and lissajous cross plot of a transient response after 
yaw-mode excitation. Rocking-mode roll-rate interference; no damping. 
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When the rolling motion of the rooking mode interferes with the response to yaw- 
mode excitation and light damping is present , as in the time histories in figure 16(a) , 
a cross plot of the roll and yaw time histories results in an oscillatory curve like that 
in figure 16 (b) . Although the curve is oscillatory , the bounding lines of the left and 
right side are parallel and constitute the bounds of the lissajous figure with damping 
present. As in figures 14(b) and 15(b) , the slope of the bounding lines is the roll-to- 
yaw ratio of the yaw mode for the tilt angle of the spring plane , 8 

sp 


p, deg/see o j 
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(a) Time histories of transient responses . 
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-.2 -.1 0 .1 .2 
IPI, deg/sec 


(b) Lissajous cross plot; |-2-| = 0.284. 

I ^ 'v]/ 

Figure 16. Time histories and lissajous cross plot of a transient response after 
yaw-mode excitation, with damping and interference from rocking mode. 
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Practical Considerations in Setting Up the Experiment 
and Exciting the Yaw Mode 

The experimental installation must be set up so that it provides rigidity in the 
suspension sling and spring attachment provisions , minimizes spring sag , minimizes 
damping, and provides adequately sensitive onboard response sensors. The method 
used to excite the yaw mode should minimize the excitation of the sway and rocking 
modes . 

The suspension sling should constitute , in effect , a rigid system between the 
vehicle and the point of suspension . The point of suspension at the end of the sus- 
pension line should constitute a single-point pivot . Single-point pivot action is 
essential to the determination of the vertical center-of- gravity position by the single- 
point suspension method . 

The pylons or outrigger assembly used to attach the springs to the vehicle 
must be rigid relative to the reference axis system . Elasticity in the assembly rela- 
tive to the reference axes results in secondary spring action in series with the 
regular springs , which makes a softer effective spring than intended. Figure 17 
shows the suspension configuration used for the HL-10 lifting body vehicle, and 



Center of gravity of — 

(a) Vehicle with test -setup 

components Inciuded 

(b) Vehicle with test-setup 

components removed 

0 Suspension beam 
0 Lead shot ballast 
0 Outrigger x 



[^Centerline of turnbuckle pin at 
\ front spring tieback 


Ax = 0.014 

I ^Spring attachment point 


Centerline of turnbuckle pin at 
rear spring tieback 


-4.108- 


-4.108- 


Figure 17. Configuration of HL-10 lifting body vehicle for the 
determination of and e. All dimensions are in meters . 
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figure 18 shows the spring attachment provisions used to swing the vehicle . The 
rigidity of the spring attachment outrigger relative to the reference axes is evident 
in the photograph . 



Figure 18. Single-point suspension of the HL- 10 E-14687 

lifting body vehicle showing spring attachment 
provisions , including spring mounting plate with 
spherical bearing . 


To enhance the freedom of action of the springs at the attachment point during 
rolling and yawing motions , a spring mounting plate with a spherical bearing was 
fabricated (fig. 18) . The tieback from the spring to the stanchion should be as 
light as possible to minimize sag in the line of action of the springs . It should also 
be longitudinally rigid to prevent secondary spring action in series with the actual 
spring . Therefore , flexible cables should not be used . Lightweight aluminum 
tubing with swaged ends is effective . 

Since the spring and the turnbuckle are heavy , care must be taken to minimize 
sag in the line from the spring attachment point to the tieback stanchion . The spring 
should be pretensioned enough to prevent obvious sag in the line of action during 
periods of motion when the spring tension is relaxed . 

Before the springs are installed , a plumb bob should be used to find the position 
of the spring attachment points on the ground . After the springs are attached and 
pretensioned , and after each change in the inclination of the spring plane , the spring 
attachment point should remain in the original position. This prevents forward or 
rearward displacement in the suspension and change in the pitch attitude , 

The vehicle should be leveled laterally as well as longitudinally by using a 
clinometer and ballast. Figure 18 shows the lead shot ballast used to level the 
HL-10 lifting body vehicle. This weight is accounted for in figure 17. 
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The electric lines that connect the onboard sensors with the remote recorders 
should be for pertinent sensors only , and they should be suspended so that damping 
of the vehicle's motion is minimized. 

The yaw mode should be excited without inadvertently exciting the pitch mode . 
Pitch rate as well as roll and yaw rates should be recorded to check for pitch-mode 
excitation . In addition , the excitation of the sway mode should be minimized and the 
excitation of the rocking mode should be minimized or avoided . At the NASA Flight 
Research Center , the yaw mode is excited by gradually building up a sinusoidal 
force that is applied at the yaw-mode frequency as far forward as possible . The 
force is positioned vertically so that the tendency to excite the rocking mode is mini- 
mized . One technique that can be used to apply the force is to extend the left and 
right index fingers laterally and apply them at the horizontal diameter of the nose 
boom or to a vertical plate attached to the vehicle. An oscilloscope arranged to show 
the lissajous figure responses of roll and yaw rates is useful for learning the tech- 
nique . The optimum application of the exciting force and the transient responses 
appears on the oscilloscope as a sloped straight line . 


Reduction of Experimental Data To Find Yawing Moment of Inertia, 

Product of Inertia, and Inclination of Principal Axis 

Data for the HL-10 lifting body vehicle are used to illustrate the reduction of ex- 
perimental data to obtain the yawing moment of inertia, product of inertia, and 
inclination of the principal axis . The portion of the experiment that deals with the 
determination of the vertical and horizontal position of the center of gravity is 
summarized in appendix A . 


Figure 17 shows the experimental configuration of the vehicle, the center-of- 
gravity positions of the vehicle for the test-^setup and clean conditions , and the 
center-of-gravity locations of the test equipment . Figure 19 shows the experl- 

and the period of oscillation as a function of 


mentally determined variation of j^j 

I ^ l \|/ 

8 for the test setup . At the point where 
sp 


= 0, the data indicate that 

I '*4 


sp 


IS 


known within ±0.1°. It should be noted that the j-^! curve is not necessarily a 

straight line . Therefore , if only a few data points had been obtained (for the outlying 
ends of the curve) , it would have been invalid to assume the presence of a hysteresis 
effect defined by straight lines depicting the upper and lower bounds of the data . 


With CO , and 6 known for the condition where i^! = 0 , the yawing moment 

'i' sp |r|^ 

of inertia, the product of inertia, and the inclination of the principal axis were deter- 
mined as shown in appendix B. The yawing moment of inertia for the experimental 
setup was obtained first. This value was then corrected for the effects of the experi- 
mental rig and the inertial effects of the apparent additional air mass . The air mass 
. effects , which were calculated using the procedures in reference 11 , are significant , 
as indicated in appendix B . 
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(b) Variation of period with 6 


Figure 19. Determination of the inclination angle, 6 , of the spring plane 

^ i I 

and the period at which ]-2-| of the yaw mode is equal to zero in the single- 

I ^ I vj; 

point suspension test of the HL-10 lifting body vehicle. 
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The following tabulation of the experimentally determined and calculated 
inertial characteristics shows that in this instance a reasonably good correlation 
was obtained for all the characteristics except the inclination of the principal axis , 
and, consequently, the product of inertia. The vehicle weight was 23,344 newtons. 
For this test condition , the experimentally determined values are considered to be 
more accurate than the calculated values . 


Quantity 

Experimentally 

determined 

Calculated 


7441.0 

7378.0 

e, deg 

4.3 

5.2 

Ixz, kg-m2 

443.0 

534.0 


CONCLUDING REMARKS 


Current methods of experimentally determining aircraft center of gravity , 
moments of inertia , and product of inertia were reviewed . 

For pitching- and rolling-moment-of-inertia determination , with the airplane 
mounted on and pivoted about knife edges , spring configurations can be set up that 
eliminate nonlinear spring moments at large amplitudes of oscillation . 


Accurate values for the yawing moment of inertia , the product of inertia , and 
the inclination of the principal axis can be obtained from yaw -mode oscillation data 
acquired from tests that use a single-point suspension double-pendulum technique 
provided that the sway-mode effects are minimized by proper design of the suspension 
rig . Rocking-mode effects in the data can be isolated . The single-point suspension 
setup also provides an accurate determination of the location of the vertical center of 
gravity . 

Care in conducting the experiment is essential. The springs must be calibrated 
accurately and pretensioned to avoid the nonlinear region of deflection at minimum 
loading conditions during perturbations . Precautions must be taken to prevent the 
introduction of undesirable secondary spring effects due to the cables or to struc- 
tural deformations . In addition , the springs must be preloaded sufficiently to avoid 
any significant sag in the line of action of the springs . Precautions must also be 
taken to insure a minimum of extraneous damping of oscillatory motion . 


Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., September 4, 1974 
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APPENDIX A 


DETERMINATION OF THE CENTER OF GRAVITY OF THE HL-10 LIFTING 
BODY VEHICLE USING THE SINGLE -POINT SUSPENSION METHOD 


The experimental setup for determining the center of gravity of the HL-10 lifting 
body vehicle by using the single-point suspension method illustrated by figure 3(a) 
is the same as the setup used to obtain the yawing moment of inertia and the product 
of inertia (figs . 10(b) and 18) , except that springs are not used and the vehicle is 
tilted in pitch by the use of loading weights . Lead shot ballast was placed on the 
right outrigger (which was used for the spring attachment in the yaw tests) as 
shown in figure 18 to level the vehicle longitudinally and laterally to within ±0.08°. 
Tape measures were fastened to the vehicle at the reference rivet points (fig. 3(a)) 
normal to the horizontal reference axis and used in conjunction with a transit to ob- 
tain the tilt angle in pitch due to the loading weights . 

The following tables summarize the detailed calculations made to obtain (1) the 
vertical center-of-gravity position of the vehicle with experimental equipment 
included and (2) the vertical and horizontal center-of-gravity positions with the 
experimental equipment removed . The loading coordinates x' and z^ and the 
distance between the rivets, 2^, are listed in figure 3(a) . 


TABLE Al.- VERTICAL CENTER-OF-GRAVITY POSITION OF EXPERIMENTAL CONFIGURATION 

[Nose loading] 


© 


© 

© 








Change in average tape 

Loading 


Front tape, m 



Rear tape , m 


reading relative to zero 







loading 


weight, 
w, N 



Decreasing 

w 

Average 

Increasing 

w 

Decreasing 

w 

Average 

Front 

Rear 

Net 




0.0050 

0.118 

0.121 

0.1185 




0 




430.36 



0.1175 

0.156 

0.159 

0.1575 

0.1125 

0.0390 

0.1515 

652.77 



0.1755 

0.171 

0.176 

0.1735 

0.1705 

0.0550 

0.2255 

875.18 



0.2285 

0 . 194 

0 . 195 

0.1945 

0.2235 

0.0760 

0.2995 

1097.59 


0.290 

0.2865 

0.210 

0.213 

0.2115 

0.2815 

0.0930 

0.3745 

1320.00 


0.341 

0.3410 

0.227 

0.227 

0.2270 

0.3360 

0.1085 

0.4445 


© 


© 


Loading 
weight , 
w, N 


= 


0 

430.36 

652.77 

875.18 

1097.59 

1320.00 


0.0511 

0.0761 

0.1010 

0.1263 

0.1500 


© 

2.964 


w 


tan 0 


3.110 

© 


w 


tan 0 


60. 861 
40.867 
30.792 
24.264 
20.733 


© 



ts 


59.524 

39.530 

29.455 

23.287 

19.396 


© 

(a) 



© = 


©X© 

24,309 


m 


1.054 

1.062 

1.060 

1.051 

1.054 


Average = 1.056 


^Results in this column are plotted in figure 3 (b) . 
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APPENDIX A (Continued) 


TABLE A2. - VERTICAL CENTER-OF-GRAVITY POSITION OF CLEAN VEHICLE 


© 

© 

© 

© 

Item 

W, N 
(a) 

Vertical distance from sus- 
pension bolt to center of 
gravity of item , m 

Moment = 0 X 0, N-m 

Aircraft with test components 

24,309.00 

1.056 

25,670.304 

Suspension beam 

-733.95 

®0.162 

“118.900 

Lead shot ballast 

-142.34 

® 1.892 

-269.307 

Two outriggers 

-88.96 

®1.892 

-168.312 


"'a/c ^ " 23.343.75 


2] moments = 25,113.785 


^From figure 3(a). 
^From table Al. 


The vertical position of the center of gravity of the clean vehicle below the sus- 
pension bolt pivot may be found as follows: 

_ y moments 

z''= 1.076 meters (Al) 

yw 

The shift in vertical position of the center of gravity due to the removal of the 
experimental gear is defined as 


Az = z'' - zl =0.20 meter 
ts 


(A2) 


The vertical position of the center of gravity below the horizontal reference iden- 
tified by reference rivets (fig. 3(a)) can be determined as follows: 

z = z^- z^ = 1.076 - 0.884 = 0.192 meter (A3) 

r r 


TABLE A3. - HORIZONTAL CENTER-OF-GRAVITY POSITION OF CLEAN VEHICLE 


© 

© 

© 

© 

Item 

W, N 
(a) 

Horizontal distance from sus- 
pension axis to the center of 
gravity of item , m 
(a) 

Moment = @ X Q), N-m 

Aircraft with test components 

24,309.00 

0 

0 

Suspension beam 

-733.95 

“0.419 

307.525 

Lead shot ballast 

-142.34 

-0.089 

12.668 

Two outriggers 

-88,96 

-0.051 

4.537 


W^/c = 2^ = 23,343.75 


2) moments = 324.730 


^From figure 3(a), 
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APPENDIX A (Concluded) 


The shift in the horizontal position of the center of gravity due to the removal of 
experimental gear can be determined from the following equation: 

_ y moments 

Ax = — — = 0.014 meter forward of suspension line (A4) 

Zw 

The fuselage station of the center of gravity of the clean vehicle can be found as 
follows: 


Fuselage station = Reference fuselage station + ~ Ax 

= 3.632 - 0.359 - 0.014 = 3.259 meters 

V ^ / 

From figure 3(a) 
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APPENDIX B 


DETERMINATION OF THE YAWING MOMENT OF INERTIA, 
PRODUCT OF INERTIA , AND INCLINATION OF THE 
PRINCIPAL AXIS OF THE HL-10 LIFTING BODY VEHICLE 
USING THE SINGLE-POINT SUSPENSION DOUBLE-PENDULUM METHOD 


Figures 10(b) and 18 show the test setup of the HL-10 lifting body vehicle for the 
determination of the yawing moment of inertia, product of inertia, and inclination of 
the principal axis . Figure 17 shows the locations of the centers of gravity of the test- 
setup components , and the location of the vehicle center of gravity (as obtained from 
appendix A) in the test-setup and clean configurations . In addition , figure 17 lists 
the linear spring constants of each of the four springs used in the setup . 

The yaw-mode variation of the roll-to-yaw ratio and the period with the inclina- 
tion of the spring plane, as obtained from test data, is summarized in figure 19, The 
only values of inclination angle of the spring plane and period of interest are those at 

which 7 ^ = 0 . 

Klv|/ 


YAWING MOMENT OF INERTIA ABOUT THE VERTICAL REFERENCE AXIS 
THROUGH THE CENTER OF GRAVITY 


Yawing Moment of Inertia, Including Test- Setup Components, 
About the Suspension Axis 

The yawing moment of inertia for the test setup is obtained from 


(8sp)„ _ (®sp)„ 


^ 0 , 

'I/ 




where , from figure 19 for 


l£l 

I T> I 


(Bl) 



3.46® 


and 


P = 1.66 sec 
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APPENDIX B (Continued) 


and , from figure 17 , 

= (1.935)^ (7442 + 7413 + 7413 + 7544) (B2) 

= 111,626 N-m/rad 


Thus , 


I = 7784.00 kg-m^ 

^s/, 

'ts 


Yawing Moment of Inertia About Suspension Axis With Test- Setup 
Components Removed and Added Air Mass Effects Accounted For 


The yawing moment of inertia of the clean vehicle about the suspension axis is 
obtained from the following equation: 


where 





I^ \ = 7784.00 kg-m^ 

®4s 


(B3) 




= 127.40 kg (obtained from the following table): 


© 

© 

© 

© 


MM 

© 

Test-setup 

component 

W, N 

© 

9.807’ 

kg 

I„ , kg-m^ 

Center of gravity 
(a) 

+ y^ , 
s s 

=©^®x®. 

0 ec 

kg-m^ 

X , m 
s 

y » m 
*^s 

Suspension beam 
Ballast 

Left outrigger 
Right outrigger 

733.95 

142.34 

44.48 

44.48 

74.84 

14.51 

4.54 

4.54 

= 46.82 

Negligible 

Negligible 

Negligible 

-0.419 

-0.089 

-0,051 

-0.051 

1 

0.176 

2.990 

2.646 

2.646 

60.00 

43.38 

12.01 

12.01 

= 12^-40kg 


^From figure 17. 
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APPENDIX B (Continued) 


The added air mass effect on ^ J , was accounted for by using the pro- 

' 'am 

cedure in reference 11 . This was calculated to be 215.57 kg-m^ . 

Thus, 


= 7784.00 - 127.40 - 215.57 = 7441.03 kg-m^ 
^s 


Yawing Moment of Inertia of the Clean Vehicle About Its Z-Axis 


The yawing moment of inertia of the clean vehicle about its Z-axis is obtained as 
follows: 


W 


h h ' g 

s ® 


a/c / . 


(a5)> 


(B4) 


where, from appendix A, = 23,343.75 newtons and Ax = 0.014 meter . Thus, 


Iz = - -^80Y 


2^34^. '^ /-n ni/i\2 -- 


(0.014)" - 7441.03 -0.47 
7441 kg-m^ 


PRODUCT OF INERTIA RELATIVE TO THE REFERENCE XZ -AXES 
THROUGH THE CENTER OF GRAVITY 


Product of Inertia , With Test-Setup Components Included , About the 

X Z -Axes Passing Through the Suspension Line 
s s 

The product of inertia for the test setup can be obtained as follows: 

tan 6 




sp 


(B5) 


where, from equations (Bl) and (B2) , n_ \ = 7784.00 kg-m^ and, from figure 19, 

IpI VVts 

f” 6^P=3.46". 


ivl/ 
Thus , 


(*X Z ) 

^ ® ® 4s 


= 470.64 kg-m^ 
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APPENDIX B (Continued) 


Product of Inertia About X Z -Axes With Test-Setup Components Removed 

s s 

The product of inertia of the clean vehicle about the X^Z^-axes is obtained from 
the following equation: 


■XoZ„ ^ 


/ee 


(B6) 


where 



= 470.64 kg-m^ 



mx z 
s s 


’)ec 


26.60 kg-m^ (obtained from the following table) : 


© 

© 

© 

© 

© 

© 

© 

Test-setup 

component 

W, N 

© 

9.807’ 

kg 

, kg-m' 


*s"s- 

kg-m^ 

X , m 
s 


Suspension beam 
Ballast 

Left outrigger 
Right outrigger 

733.95 

142.34 

44.48 

44.48 

74.84 

14.51 

4.54 

4.54 

Negligible 

Negligible 

Negligible 

Negligible 

M 

-0.894 

0.836 

0.836 

0.836 

0.3750 

-0.0744 

-0.0426 

-0.0426 

28.06 

-1.08 

-0.19 

-0.19 

V/ly y + iTfix = 26.60 kg-m^ 

2,\ ^oh ® 7eo 


^From figure 17 . 


Thus, 


I 


X Z 


s s 


= 470.64 - 26.60 


444.04 kg-m^ 


Product of Inertia of the Clean Vehicle About Its XZ-Axes 

The product of inertia of the clean vehicle about its XZ-axes is obtained as 
follows: 


^XZ ^X Z 
s s 



(B7) 
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APPENDIX B (Concluded) 


where , from table A2 , 

= 23, 343 . 75 newtons 

Az = 0 . 020 meter 
and , from table A3 , 

AjE = 0.014 meter 
Thus , 

Ixz = 444.04 - (0.014) (0.020) = 443.4 kg-m^ 

INCLINATION OF THE PRINCIPAL AXIS 


With Iz 

443 . 4 kg-m^ 

1625.0 kg-m 
tion (25): 


and I^z’ obtained in this appendix , equal to 7441.0 kg-m^ and 
, respectively , and I^, as obtained from other sources , equal to 
^ , the inclination of the principal axis can be obtained by using equa- 


from which 


tan 2s 




0 . 1525 


e = 4.34® (below the reference X-axis) 
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